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^ ■ Abstract 

d ' It has been pointed out by Shlyakhter that data from the natural 

fission reactors which operated about two billion years ago at Oklo 
(Gabon) had the potential of providing an extremely tight bound on 
the variability of the fine-structure constant a. We revisit the deriva- 
tion of such a bound by: (i) reanalyzing a large selection of published 
rare-earth data from Oklo, (ii) critically taking into account the very 
large uncertainty of the temperature at which the reactors operated, 
and (iii) connecting in a new way (using isotope shift measurements) 
the Oklo-derived constraint on a possible shift of thermal neutron- 
capture resonances with a bound on the time variation of a. Our final 
(95% C.L.) results are: -0.9 x 10~ 7 < (a° kl ° - a QOW )/a < 1.2 x 1(T 7 
and -6.7 x lO^yr" 1 < d avcraged /a < 5.0 x lO^yr" 1 . 
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1 Introduction 



Since Dirac ]T] first suggested it as a possibility, the time variation of the 
fundamental constants has remained a subject of fascination which motivated 
numerous theoretical and experimental researches. For general discussions 
and references to the literature see, e.g., [@, ^ |], |E|. Superstring theories 
have renewed the motivation for a variation of the "constants" by suggesting 
that most of the dimensionless coupling constants of physics, such as the 
fine structure constant a = 1/137.0359895(61), are functions of the vacuum 
expectation values of some scalar fields (see, e.g., ||). Recently, a mechanism 
for fixing the vacuum expectation values of such massless stringy scalar fields 
(dilaton or moduli) has been proposed 0. This mechanism predicts that 
the time variation of the coupling constants, at the present cosmological 
epoch, should be much smaller than the Hubble time scale, but maybe not 
unmeasurably so. In this model, the time variations of all the coupling 
constants are correlated, and the ones of most observational significance are 
the fine structure constant a and the gravitational coupling constant G. In 
the present paper, we revisit the current best bounds on the variation of a. 

One of the early ideas for setting a bound on the variation of a was to 
consider the fine-structure splittings in astronomical spectra |3J. With this 
method, Bahcall and Schmidt concluded that a had varied by at most a 
fraction 3 x 10 -3 of itself during the last 2 x 10 9 years. A recent update of this 
method has given the result Aa/a = (0.2 ±0.7) x 10~ 4 at redshifts 2.8 < z < 
3.1, i.e. the bound \a/a\ < 1.6 x 10~ 14 y r_1 (2cr level) on the time derivative 
of a averaged over the last ~ 10 10 yr [III. See also Ref. [TT| which obtains 



—4.6 x 10~ 14 yr~ 4 < a/a < 4.2 x 10~ 14 yr~ 4 from fine-structure splittings, 
and, denoting x = a 2 g p (m e /m p ), —2.2 x 10 _15 yr _1 < x/x < 4.2 x 10 _15 yr _1 
by comparing redshifts obtained from hyperfine (21cm) and optical data. 

One of us obtained the upper limit \a/a\ < 5 x 10~ 15 y r 1 from an analysis 
of the abundance ratios of Rhenium and Osmium isotopes in iron meteorites 
and molybdenite ores 0. The most recent direct laboratory test of the 
variation of a has obtained \a/a\ < 3.7 x 10 -14 yr~ 4 by comparing hyperfine 
transitions in Hydrogen and Mercury atoms For more references on the 
variation of constants see 0, ^ 

On the other hand, the much more stringent bound \a/a\ < 10 _17 y r_1 
has been claimed by Shlyakhter |1^, |15[] to be derivable (at the three stan- 
dard deviations level) from an analysis of data from the Oklo phenomenon. 
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The Oklo phenomenon denotes a natural fission reactor (moderated by wa- 
ter) that operated about two billion years ago in the ore body of the Oklo 
uranium mine in Gabon, West Africa. This phenomenon was discovered by 
the French Commissariat a l'Energie Atomique (CEA) in 1972. The results 
of a thorough, multi-disciplinary investigation of this phenomenon have been 



presented in two conference proceedings fl6 , [T?| . See also Jl8| and [19| for 
summaries of the first phase of investigation. 

In view of the importance of Shlyakhter's claim, of the lack of publication 
of a detailed analysis Q, and of our dissatisfaction with some important as- 
pects of the analysis presented in two preprints |14|, |T5j| , we decided to revisit 



the Oklo bound on a. The main conclusions of our work are the following: (i) 
we confirm the basic claim of Shlyakhter that the Oklo data is an extremely 
sensitive probe of the time variation of a; (ii) after taking into account var- 
ious sources of uncertainty (notably temperature effects) in the analysis of 
data, and connecting in a improved way the raw results of this analysis to 
a possible variation of a, we derive what we think is a secure (95% C.L.) 
bound on the change of a: 



-o.9xi(r 7 < <i.2xi(r 7 . (i) 



„Oklo „,iiow 

a 

In terms of an averaged rate of variation, this reads 

.-17,-1 „ a ^ r n in-17, -J 



- 6.7 x Kn 17 yr < - < 5.0 x fCr^yr" 1 • (2) 



a 



Extracting the neutron capture cross sec- 
tion of Samarium 149 from Oklo data 



The proof of the past existence of a spontaneous chain reaction in the Oklo 
ore consists essentially of: (i) a substantial depletion of the Uranium isotopic 
ratio tj235^tj238 reS p ec t to the current standard value in terrestrial 
samples; and (ii) a correlated peculiar distribution of some rare-earth iso- 
topes. The rare-earth isotopes are abundantly produced in the fission of 

1 The very brief account published in Nature |TJ| omits most of the analysis that is 
presented in the two preprints jl4], p!B| . 

2 For typographical convenience, we indicate atomic mass numbers as right, rather than 
left, superscripts. 
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U 235 and the observed isotopic distribution is beautifully consistent with cal- 
culations of the effect of a strong neutron flux on the fission yields of U 235 
(see e.g. GO, |STL 1231, E3L |2H|). In particular, the strong neutron absorbers 



Sm 149 , Eu 151 , Gd 155 and Gd 157 are found in very small quantities in the cen- 
tral regions of the Oklo reactors (see, e.g., Fig. 2 of p3l). These isotopes 



were evidently burned up by the large neutron fluence produced by the fis- 
sion process. Following Shlyakhter's suggestion [14, 15fl , we concentrate on 



the determination of the neutron capture cross section of Sm 149 : Sm 149 (n, 7) 
Sm 150 . 

The evolution of the concentrations of the various Samarium isotopes 
(sharing the common atomic number Z = 62) in the Oklo ore is especially 
simple to describe because of the absence of a stable chemical element with 
atomic number Z = 61. The most stable Promethium nuclide is PrnJ 45 with 
a half-life of 17.7 years. If there had existed, before the reaction started, 
some natural concentration of Promethium it could, via neutron absorption 
and subsequent (3~ decay, have generated some Samarium. In absence of 
this, the final values of the Sm concentrations are determined by: (i) their 
initial concentrations, before the reaction; (ii) the yields from the fissions; 
and (iii) the effect of neutron captures. 

Following Refs. (16], [17|], one characterizes the neutron absorbing power 
of an isotope by the effective cross section 

_ fa(E)vn E dE 
v J uecLE 

where v is the (relative) velocity of incident neutrons, n E dE the energy dis- 
tribution of the neutrons, and vq the fiducial (thermal) velocity vq = 2200m/s 
corresponding to a kinetic energy Eq = 0.0253eV. The advantage of the defi- 
nition (U) is that, in the case of a u l/v absorber", a(E) = C/v, the effective 
cross section equals a = C/vq = <j(Eq) independently of the neutron spec- 
trum. [To a good appoximation, this is the case for the thermal fission cross 
section of U 235 ] . On the other hand, in the case of nuclides exhibiting reso- 
nances in the thermal region (these are the strong absorbers, Sm 149 , Eu 151 , 
Gd 155 , Gd 157 ), the value of the effective cross section (|3|) is very sensitive to 
the neutron spectrum, especially to its thermal part Q 

3 The spectrum of moderated neutrons in a fission reactor consists of a Maxwell- 
Boltzmann thermal distribution up to energies of order a few times kT, followed by a 
tail nsdE oc dE/(vE) due to neutrons still in the process of moderation 
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Associated to the introduction of the effective cross section (|3|), one de- 
fines an effective neutron flux = uvq with n = J hecIE, and an effective 
infinitesimal fluence (integrated flux) 

dr = cf)dt = nvodt . (4) 

With this notation, the general equation describing the evolution of the total 
number N A of nuclides of mass number A (for some fixed atomic number Z) 
in some sample reads 

-j-^ = VAN 5 a f5 + a A ~iN A -i - <?aN a ■ (5) 
ar 

Here, y A denotes the yield of the element A in the fission of U 235 , N§ and 
ays are short hands for A^ 235 and the (effective) fission cross section of U 235 , 
and the last two terms describe the effects of neutron captures within isotopes 
of some chemical element Z. For simplicity, we drop the carets over the cross 
sections. The evolution equation (|) neglects any contribution oc N A -i(Z') 
coming from the (3~ decay of the neighbouring chemical element Z' = Z — 1 
after absorption of a neutron. As we said above, this approximation applies 
well to the Samarium case. Eq. @ neglects also the yields due to the frac- 
tionally small number of fissions of U 238 and Pu 239 . [See, e.g., Ref. [2~3"| which 



estimates that, in a particular sample, 2.5% and 3% of the fissions were due 
to U 238 and Pu 239 , respectively]. 

Samples in the cores of the various Oklo reactors were exposed to a total 
effective fluence r — J dr = J nv^dt of the order of 10 21 neutron/cm 2 = 1 
inverse kilobarn. This means, roughly speaking, that processes with effective 
cross sections comparable or larger than 1 kb have led to a significant number 
of reactions, while processes with a <C lkb had a negligible effect. The former 
category includes the fission of U 235 (cr/ 5 ~ 0.6kb), the capture of neutrons by 
Nd ((7i43 ~ 0.3kb) and by the strong absorbers (such asSm 149 ;ai4 9 ^70kb), 
while the latter category includes neutron captures by weak absorbers such 
as Sm 14 and Sm 148 with cross sections of only a few barns. 

This allows one to neglect o"i 44 and a U8 (for Z = 62) in Eq. (|). Further 
simplification comes from the fact that the stable isotopes 144, 146 and 148 
of Neodymium prevent the formation of the long-livedP| Sm 144 , Sm 146 and 

4 The half-life of Sm 146 is 1.03 x 10 8 yr and therefore long with respect to the duration 
of the Oklo phenomenon. 
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Sm 148 as end points of (/3 decay) fission chains. The stable isotopesQ obey 
the simple evolution equations 



tZJV, 



l-M 



dr 
dN 147 

dr 
dNus 

dr 
dN 149 

dr 
dN 5 

dr 



0, (6) 

2/147^5^/5 - ^147^147 , (7) 

^147^147 , (8) 

2/149^50/5 - 0-149^149 , (9) 

-N s o$ . (10) 



To close the system, we have followed [^2j and |23j in describing the burn 
up of U 235 by means of a modified absorption cross section a£ = a 5 (l — C), 
where 05 is the normal absorption cross section (fission plus capture) and C 
is a conversion factor representing the formation of U 235 from the decay of 
Pu 239 formed by neutron capture in U 238 . 

In the approximation where the (effective) cross sections, and the conver- 
sion factor C, are constant, the system fl6HlO|) is easily solved and gives 

N b (t) = N 5 (p)e-°i T , (11) 
iV 144 (r) = iV 144 (0), (12) 

N U7 {t) + iV 148 (r) = iV 147 (0) + iV 148 (0) + ywaf 5 N 5 (0) 1 ~ ^ - , (13) 

f>~ a t T g — fT 149 T " 

N U9 (t) = jV 149 (0)e- CTl49T + 2/i49a /5 iV 5 (0) — j^- .(14) 

0"l49 — 0"5 

Eq. (0) shows that the quantity of Sm 144 measured in a sample now is equal 
to the quantity of natural Sm 144 present in the sample before the nuclear 
reactions. Assuming that the natural Samarium present in the sample at the 
beginning had the normal isotopic ratios (say n 144 = 3.1%, n 147 = 15.0%, 
nus = 11.3%, ni 49 = 13.8%, etc... [PEfl), we can use flT2"|) to correct flT3|) 
for the initial concentrations in Sm 147 and Sm 148 . The effect of A^i 4 g(0) 

5 "Stable" means, in this context, a half life much larger than the age of the Oklo 
phenomenon. E.g. the half-life of Sm 147 is 1.06 x 10 11 yr^> 2 x 10 9 yr. 
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in Eq. (|TJ]) is totally negligible (as is the last term) because the exponent 
(T149T 3> 1. We then derive the intermediate result 



, in r) + iV 148 (r) - BM2±m^ NlM{r) y ^ ^ _ x 



>149 - CT5) • ( 15 ) 



One can finally obtain an expression for cr 149 in terms of "measured" 
quantities by connecting a\r to the observed ratio between the numbers of 
U 235 and U 238 atoms in the Oklo sample. If, following we define 



w 



0.00725 

we 



(16) 



where 0.00725 is the usual tj235^tj238 ra ^j j n na tural Uranium now, it is easy 
to verify that w = e a $ T . Finally, we get Q 



C"l49 



1 

T 



In w + y 



'148 



nN- 



111 



w 



N- 



119 



(17) 



where (using Refs. [^7| and [J 

2/149 



1.080384 



V 



11 



0.478 



y 147 2.261681 
rii4 7 + ni4 8 15.0 + 11.3 



71144 



3.1 



8.48 . 



(18) 
(19) 



The quantities Na in Eq. ([i~7|) denote the present values of the isotopic 
concentrations, or, equivalently, the present values of the isotopic ratios. 
The isotopic ratios of Samarium have been measured in many Oklo sam- 
201, E31 I 



pies 



Note that the quantity which is, at this stage, 
directly obtainable from observations is the dimensionless product (J\i\qt = 
J dt J <Ti49(E)vriEdE. 

Similarly, by considering the fission yields of Neodymium and the neutron- 
capture reaction Nd 143 — >Nd 144 (using, e.g., the Nd 142 content to subtract the 
contribution from the natural concentrations present before the reaction), 

2 , [l4 , ll| apart from the facts 
l)/lnu; instead of its inverse 

in Eq. (p7|)). The fractionally small first contribution on the right-hand side of ( |l7| ) is 
neglected in the above references. 



6 Eq. ( J17| ) is equivalent to equations appearing in Refs 
that Shlyakhter's equations contain misprints (e.g. (w 
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several authors |20], ^T], ^] have shown how to obtain the dimensionless 
product <7i43T (where 0143 = <j( ni7 )(Nd 143 )) in terms of quantities observed in 
Oklo samples. Combining these two results, we see that the value two billion 
years ago of the ratio (Juqt / ai 43 r = <7i 4 g/cri 4 3 can be computed in terms of 
present Oklo data. 

Although it would be conceptually clearer to deal only with the dimen- 
sionless ratio cr 149 /cTi43, we shall follow previous usage in working with the 
dimensionful quantity o"i 4 g obtained by inserting in Eq. (|T7|) the value of 
the (effective) ffuence r deduced from Neodymium data by previous authors. 
This procedure is justified by the fact that the effective cross- section 17143, de- 
fined by Eq. (||), depends very little on the neutron spectrum because a\^{E) 
follows the 1/v law over most of the range of interest. Therefore the lack of 
knowledge of the temperature of the moderated neutrons is of no importance 
(contrary to the case of cr^g) and the effect of epithermal neutrons is also 
very smallQ In other words, the extraction of CT143T from Oklo data is approx- 
imately done by assuming a fixed, fiducial value for (X143, say (J143 ~ 325b, so 
that the use of Eq. ( JT7| ) for computing a dimensionful <Ti 4 g is approximately 
equivalent to computing the dimensionless quantity 325cri 4 g/<7i43. 

The detailed isotopic analysis of Oklo data [JR], |T7]] has shown that, gen- 
erally speaking, the ore composition has changed very little since the end of 
the nuclear reactions. This is established by studying the correlation between 
the ffuence r and the Uranium isotopic ratio N^/Ng, and by showing that 
it can be explained by neutronics considerations (see, e.g. [p9p . However, 



in some cases there is evidence for a partial reshuffling of chemical elements 
after the end of the reactions. We have examined these results and selected 
16 samples as especially suitable for extracting a reliable value of <Ti4g. These 
samples are all core samples with high Uranium content, large depletions of 
U 235 , and high fluences, r^0.7 x 10 21 n/cm 2 . In all cases, the natural ele- 
ment correction in Eq. ( |17D is small, the observed Samarium having been 
produced almost entirely by fission. The very small content of Sm 149 (and, 
when data are available, of other strong absorbers such as Gd 155 and Gd 157 ) 
is also a confirmation of the absence of chemical reshuffling after the reaction. 
The data we took come from [28], 22], |29|, |24|], and |31|0. The result of 



'In the analysis of Oklo data, it has been customary to parametrize the contribution 
of epithermal neutrons to the spectrum by a parameter called r. This parameter is found 
to be small, r ~ 0.15, and its effect on 0143 is only a few percent f3(i|| . 

8 The sample SC521472 taken from this last reference was exposed to a smaller fluence 
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calculating aug from these data is exhibited in Table I. 

The large scatter of the values exhibited in Table I is compatible with 
the strong temperature dependence of a U9 (see below). The only exception 
is the 36kb obtained for the sample SC39-1387. This value is a clear outlier 
which, most plausibly, has been contaminated in some way. Excluding this 
result, the other 15 results are all contained in the range 

57 kb < <7 149 < 93 kb . (20) 

We think that it is conservative to consider the full range (^) as a "2cr" (or 
95% C.L.) interval for a^g. [For clarity, we reestablish the caret meaning 
that we are dealing with an effective cross section, Eq. (0).] Actually, in 
view of what is known from Oklo, it is very plausible that the range (pC|) is 
to be attributed to a mixture of temperature effects and a small amount of 
post-reaction chemical reshuffling. For our purpose, we will use the full range 
( |2"U| ) to define a conservative bound on the variation of a. For completeness, 
and in view of the special use we make of the interval (p0|) we give in Table II 
the complete set of data allowing one (using (|T7D) to compute a^g for the 
samples giving the extreme values (|20|). 

Let us note that the values we obtain for a^g are different from the result 
claimed by Shlyakhter [14], [15[], namely a^g = (55 ± 8)kb. As he did not 



mention the data he used, we could not trace the origin of this difference. 
We note that most of the values in Table I are compatible with thermal 
effects (aug increases from ~ 70kb to ~ 99kb when the temperature varies 
from 20°C to ~ 400°C, and then decreases for higher temperatures^) . 

3 Bounding a possible shift of the lowest res- 
onance in the capture cross section of Samar- 
ium 149 



Following Shlyakhter's suggestion [O, |T3, ITS], we shall translate the range of 



'Oklo" values of a U g, Eq. (£(])> i n t° a bound on the possible shift, between 



than the others. It was included because this sample has been used to estimate the 
temperature of the neutrons. 

9 The fact that we did not find values between 93 and 99kb is probably explained by 
some post-reaction remobilization. Anyway, the limit we shall derive on a depends only 
on the lower bound on (T149. 
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the time of the Oklo phenomenon and now, of the lowest resonance in the 
monoenergetic cross section 0149 (£7). The large values of the thermal capture 
cross sections of Sm 149 , Gd 155 and Gd 157 are due to the existence of resonances 
in the thermal region. In presence of such a resonance, the monoenergetic 
capture cross section is well described, in the thermal region, by the Breit- 
Wigner formula 

, s h 2 TJE)T~, , . 

< 21) 

Here p is the momentum of the neutron, E = p 2 /2m n its kinetic energy, 
g = (21' + l)(2s+ 1)~ 1 (2I + 1) -1 a statistical factor depending upon the spins 
of the compound nucleus I', of the incident neutron s = |, and of the target 
nucleus /, T n (E) is a neutron partial width, T 7 a radiative partial width, and 
r the total width. The neutron partial width T n (E) varies approximately as 

£1/2 

F n (E) = ^p , (22) 

where 7^ is a "reduced partial width" (see, e.g., [|32|). With sufficient ap- 
proximation, the total width is given by 

r ~ r 7 + r n (E r ) . (23) 

As we shall see explicitly below, the position of the resonance E r (with 
respect to the threshold defined by zero-kinetic-energy incident neutrons) is 
extremely sensitive to the value of the fine-structure constant. By contrast, 
the other quantities entering the Breit-Wigner formula, 7^, T 7 , have only a 
mild (polynomial) dependence on a. Therefore the sensitivity to a of the 
effective cross sections a, Eq. (|), measured with Oklo data is totally domi- 
nated, for strong absorbers, by the dependence of u(E) upon the position of 
the lowest lying resonance E r . As we said above, we should, more rigourously, 
work with dimensionless ratios such as aug/ai^. However, the mild absorber 
Nd 143 has no resonances in the thermal region. Therefore the a-sensitivity 
of the ratio (T149/CT143 is completely dominated by the a-sensitivity of &ug 
inherited from the dependence on E^ A9 (a). 

The main problem is to use the range (|TJ) of values of (J149 to put a limit 
on a possible shift of the lowest lying resonance in Sm 149 , 

A = £ r 149 (° kl °) - £ r 149 ( now ) . (24) 
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Previous attempts J[3j, [b|, [L^, |J at relating Oklo-deduced ranges of values of 
(Tug to A are unsatisfactory because they did not properly take into account 
the very large uncertainty in a U9 due to poor knowledge of the neutron tem- 
perature in the Oklo reactors. The original analysis of Shlyakhter assumed 
a temperature T ~ 20°C (which is much too low), and the analysis of [|J 
took T ~ 1 000 K, i.e. T ~ 725°C (which is possible, but on the high side) 
and assumed that one could work linearly in the fractional shift A/E r . We 
think that one should neither fix the neutron temperature T (which could 
have varied over a wide range) , nor work linearly in A/£ r (which could have 
been larger than unity). 

Several studies, using independent data, have tried to constrain the value 
of the temperature in the Oklo reactors W7\. Mineral phase assemblages 



observed within a few meters of the Oklo reactor zones 2 and 5 indicate a 
minimum temperature in these regions of about 400°C, while relict textures 
in the reactor zone rock suggest that temperatures T ~ 650 — 700° C may 
have been reached within the reactors |33fl . A study of fluid inclusions and 
petrography of the sandstones suggest pressures p ~ 800 — 1000 barQ and 
temperatures ranging between 180°C and at least 600°C On the other hand, 
the temperature of the water-moderated neutrons during the fission reactions 
has been evaluated by a study of the Lu 176 /Lu 175 and Gd 156 /Gd 155 isotope 
ratios in several samples. The values obtained range between 250±40°C and 



149 



450±20°C depending upon the sample and the isotope ratio considered [[H 
It is to be noted that the concentration of strong absorbers such as Sm 
or Gd 155 (which are burned very efficiently) is determined by the values of 
the effective cross sections a U9 or <r 155 at the end of the fission phenomenon. 
Therefore, we cannot exclude that the temperature to be used in evaluating 
<Ti49 or (Ti55 be on the low side of the allowed range. 

Summarizing, we consider that the temperature to be used to determine 
(5"i49 or <5"i55 could be in the full range 

180°C < T<700°C . (25) 



10 It is thought that the Oklo phenomenon took place while the Uranium deposits were 
buried ~ 4 km deep p3|. 
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Inserting a Maxwell-Boltzmann spectrum R 

n f dE = jJk* e ^ El,2dE • (26) 

in the definition (|3]), with <j(E) of the Breit-Wigner form (pl~D, we find that 
the dependence of the effective cross section of a strong absorber on resonance 
shift A, Eq. (|24T), and temperature is given by 



a(A,T)= ,^ /2 a (l + ?/ 2 ) / (27) 



(vrfcT) 3 / 2 uv yu/ 7o l + 2/ 2 (£,A) 

Here <To denotes the "thermal" radiative cross section (as observed now), 
i.e. the monoenergetic <7( ni7 )(£b) evaluated at E = 0.0253eV, yo denotes 
2(£ -£ r now ) /r, and y(E, A) = 2(E-E™ W - A) /r. The £ 1/2 in the integrand 
comes from combining several different factors: a factor E~ 1 / 2 coming from 
T n (E)/p 2 (1/v law), a factor E 1 ! 2 coming from the factor v in Eq. (|3]), and 
the factor E 1 ! 2 in the Maxwell spectrum (|26|). 

In the case of Sm 149 the numerical values needed to evaluate (^7]) are 
(from |§): E r now = 0.0973eV, a = 40.14kb, and T ~ 0.061eV. [The 



latter being estimated from T ~ T 7 + r n (£? r ) with T 7 = 60.5 x 10~ 3 eV, 
2gT n (E r ) = 0.6 x 10 _3 eV, with 2g = 9/8 corresponding to I' = 4 and 
/ = 7/2.] The dependence of o"i 4 g upon the resonance shift A is shown in 
Fig. 1 for several temperatures spanning the range (|2"5|). On the same Figure, 
we have indicated the conservative range of values of 0149, Eq. (p0|). 

The limits on A shown in Fig. 1 are — 0.12eV< A < 0.08eV. The lower 
limit depends on the minimal allowed temperature. Given the temperature 
estimates quoted above, we consider 180°C as a firm minimal temperature 
and therefore —0.12 eV as a firm lower bound. As the upper limit 0.08eV 
depends on the maximum allowed temperature which is more uncertain, we 
have also explored temperatures higher than 700° C. We found that when 
A = 0.09eV a never exceeds 57kb even if T is allowed to take values much 
larger than 700°C [cr(0.09, T) reaches a maximum < 57kb somewhere around 
T ~ 1000°C.l Therefore to be conservative, we take A < 0.09eV as firm 



11 We do not consider the effect of epithermal neutrons, which introduce only a rather 
small fractional correction (spectrum index of order r ~ 0.15 17|). This correction is 
negligible compared to the wide range we consider. 
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upper limit. We conclude that the Oklo Samarium data constrain a possible 
resonance shift to be in the range 



0.12 eV < A < 0.09 eV . (28) 



For comparison, let us mention that Refs |13], 14, pl| estimate a 2a range 



|A| < 0.02eV from the Samarium data alone, and a 3<r range |A| < 0.05eV 
from combining Samarium and Europium data. 

We tried to make use of Oklo Gadolinium data to restrict further the 
Samarium-derived range fl2"8"|). A priori, one could think of making use of 
both Gd 155 and Gd 157 which are strong absorbers of neutrons. In fact, Gd 157 
is such a strong absorber that its final concentration ^157 oc 2/157/0-157 (gen- 
eralizing Eq. ([14]) ) is too small to be measured reliably. The case of Gd 155 
is more favorable, its effective cross section being comparable to that of 
Sm 149 . However, its fission yield is much smaller (1/155 = 0.032% instead of 
?/i49 = 1.08%). The absolute concentration of all isotopes of Gadolinium is 
about ten times smaller than that of Samarium (see, e.g., p3[). This im- 
plies that Gadolinium data are much more prone to various contaminations 
(natural element contamination due to a post-reaction remobilization, and 
uncertainties in the isotopic analysis measurements). To make a meaningful 
analysis of Gadolinium data, one should probably restrict oneself to samples 
that were exposed to rather mild fluences. Such samples are SC361901 and 
SC521472 which have been studied in detail in The effective cross sec- 
tions a(Gd 155 ) obtained in the latter reference are 0-155 = (42.0 ± 0.5)kb in 
SC 361901, and 0155 = (32.5±0.5)kb in SC 521472. These values are compat- 
ible with the present values of 0-155 if the temperatures in these samples were 
T 36 i9oi — 380°C and T521472 — 450°C. [Actually, these temperatures disagree 
with the Lutetium-derived ones: T 3 L 6 U 1901 ~ 250°C and T 5 L 2 U 1472 ~ 280°C. This 
difference is probably to be explained by a moderate amount of contamina- 



tion of natural Lutetium after the reaction |JT|.] However, we could not use 
these data to derive more stringent limits on a because the Gd 155 resonance 
turns out to be less than half as sensitive as Samarium to changes in a (see 
next Section). 



13 



Translating possible resonance shifts into a 
bound on the variation of the fine-structure 
constant. 



Let us finally translate the allowed range ( |28| ) into a bound on a possible 
difference between the value of a during the Oklo phenomenon and its value 
now. The treatments given in previous analyses are unsatisfactory. The orig- 
inal analysis of Shlyakhter |14], |U| rested on a coarse representation of the 



nucleus as a square potential well, together with dubious assumptions about 
nuclear compressibility, while the analysis of Ref. [[| used an ill-motivated 
finite-temperature description of the excited state of the compound nucleus. 

The observed neutron-resonance energy E^ ow = 0.0973 eV, for the radia- 
tive capture of neutrons by Sm^ 9 , corresponds to the existence of a particular 
excited quantum state of Sm^ - More precisely, if we write the total mass- 
energy of the relevant excited state of Sm^ as Q -E* 50 = 62m p + 88m n + Ei 
(with Ei < 0), and the total mass-energy of the ground state of Snig| 9 as 
Eu9 = 62m p + 87m n + E 2 (with E 2 < 0), we have 

E r = El 50 — E U9 — m n = E\ — E 2 . (29) 
Both Ei and E 2 are eigenvalues of the Hamiltonian 

H = H n + H c , (30) 

where H c is the Coulomb energy 

^ c = e 2 E^ 1 , (31) 

summed over the pairs of protons in the nucleus, and H n is, to a good ac- 
curacy, independent of e 2 . We neglect small effects such as the magnetic- 
moment interactions or the QED corrections to the masses, such as m n and 
m p , entering the nuclear Hamiltonian H n . 

Now let e 2 vary while H n remains fixed. Then for any eigenstate of H 
with eigenvalue E, 

e 2 dE/de 2 = (H c ) , (32) 

12 We set c = 1. 
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and therefore 

e 2 dE r /de 2 = {H c ) 1 - (H c ) 2 . (33) 

The Coulomb energies on the right of (|33|) are not directly measurable. The 

quantities that can be directly measured by optical spectroscopy are 



the mean-square radii (r 2 ) of the charge-distributions of the protons in the 
various isotopes of Samarium. Let us recall that "isotope shifts" in heavy 
atoms are related to the effect of the finite extension of the nucleus on electron 
energies. A first-order perturbation analysis of the latter effect (see, e.g., f37|) 
yields AE = (27r/3)^(0)Ze 2 (r 2 ) where ip e (r) is an (s-state) electron wave 
function, and where (r 2 ) = Z^ 1 J pr 2 dv with p denoting the proton charge 
distribution in the nucleus. 

To connect the expectation-values in (|33"D with the mean-square radii, we 
use the semi-classical approximation 

(H c )i = ^e 2 Jv iPi dv, i = 1,2 , (34) 
where pi is the density of protons in the nuclear state i, normalized to 

Pi dv = Z = 62 , (35) 



for Samarium, and Vj, is the electrostatic potential generated by p^. From 
51 and (131 



dE r /de 2 = ^J {Vipi-V 2 p 2 )dv = - J 5V5pdv + J V x 5pdv < J V x 8pdv , 

(36) 



with 



5p = Pl -p 2 , 5V = V l -V 2 . (37) 



The term that is dropped in ([$6|) is negative because it is minus an electro- 
static self-energy. The integrand on the right side of (|36]) is a small difference 
5p multiplied by the smooth potential V\. With an error of second order in 
small quantities, we may approximate V\ by the classical potential of a uni- 
formly charged sphere with radius R\, 

I"-?/? 2 _ r 2 l 

y >« = 2 ^sr • (38) 



15 



Then 



where 



becomes 



dE r 
~d& 



< 



2R\ 



6 l2 (r 2 ) 



Si 2 (r 2 ) = ^ J r 2 5pdv 



(39) 



(40) 



is the difference in mean-square charge-radius between the states 1 and 2. 
Let the label 3 denote the ground-state of Sm 150 . Then 



6 l2 (r 2 ) = 5 13 (r 2 ) + 6 32 (r 2 ) 



(41) 



The difference S\ 3 (r 2 ) cannot be calculated because we do not know the shape 
of the excited state of Sm 150 . But it seems safe to assume that the proton 
charge distribution will not be more tightly concentrated in the excited state 
than in the ground state. That is to say, 

5 13 (r 2 ) > . (42) 

An inequality stronger than fl42|) could be deduced from more dubious 
assumptions about nuclear compressibility, but a stronger inequality is not 
needed. From (|39[) and (T2) we have 



dE r 
~dl? 



< 



2R\ 



Ss2(r 2 ) 



(43) 



and this is sufficient for our purposes. The experimental isotope-shift mea- 
surements reported by |38j give directly 5 3 a{t 2 ) = 0.303 ± 0.016 fm 2 and 
^24(^ 2 ) = 0.092 ± 0.005 fm 2 , where the label 4 denotes the ground-state of 
Sm 148 . Taking the difference gives 

5 32 ( r 2 ) = 0.211 ± 0.017 fm 2 (3a error) . (44) 

For the radius of the Sm 150 nucleus to insert in (15), we use equations (50) 
and (51) on page 568 of |36j, which give 

Rx = 8.11 fm. (45) 

From (B3), ( 



and (f4"5D, we find 
dR 



a- 



da 



< -(1.09 ±0.09) MeV . 



(46) 
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The estimate (^6]), obtained here directly from measurements of the small 
charge-radius difference (]44]) between Sm 150 and Sm 149 , agrees with the re- 
sult obtained by differencing the phenomenological Bethe-Weizsacker formula 
(droplet model). The latter formula estimates the nuclear Coulomb energy 
as (H c ) = 0.717 Z(Z - l)^ 1 / 3 MeV |3|. Taking the difference between 
Smgl and Snig 49 (and arguing as above that excited states are less charge 
concentrated) yields the inequality adE r /da < —1.14 MeV, which is com- 
patible with the result ((46"|). By contrast, the Bethe-Weizsacker formula 
overestimates by about a factor two the a-sensitivity of the resonance en- 
ergy E™ w = 0.0268 eV, for the radiative capture of neutrons by Gd^ 5 . This 
follows (using the same method as above) from the fact that isotope-shift 
measurements reported in [H0[ yield 



(r 2 )i56 - (r 2 )i55 = 0.097 ± 0.005 fm 2 , (47) 
which is less than half the Samarium difference (£1). As a consequence the 



a-sensitivity parameter \adE r /da\ of the Gd 155 resonance is less than half 
that of the Sm 1 9 resonance. 

We are now in position to convert the bound fl28|) obtained above from our 
analysis of Oklo data into a bound on the variation of a. To be conservative 
we use the worst 3a limit on the a-sensitivity of E r obtainable from (|46|). 
namely 

<ltr > (1.09 - 0.09) MeV = 1.0 MeV (48) 



da 



Combining (|48| ) with the bound (|^) on the shift A = E r (a° kl °) 
-\adE r /da\(a° kl ° - a now )/a yields our final result 



0.9 x 10~ 7 < — < 1.2 x 10~ 7 , (49) 



..Oklo 

a: 

which we consider 95% C.L. limit. 

Though we have been very conservative in our analysis, our result (f!9|) 
confirms the main claim of Refs. |T3| , |T4| , [T5] , namely that Oklo rare-earth data 
are extremely sensitive probes of a possible variation of a: the surprisingly 
good ~ 10~ 7 bound comes mainly from the 10 7 amplification factor between 
the MeV level in adE r /da (which is physically clearly understood from the 
Bethe-Weizsacker formula) and the 0.1 eV level of the value of the Samarium 
resonance (with respect to the threshold). 
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Let us note also the consistency of the approximations we made: a change 
Sen/ 'a ~ 1CT 7 has a totally negligible effect in all the quantities (such as 7„ 
or T) depending at most polynomially on a, and its effect on the dimen- 
sionless ratio cr U g/ai4 3 is dominated by the Breit-Wigner denominator of 
<t 149 , i.e. by the change <9£y 49 /da 8a. [The contribution coming from the 
shift of neutron capture resonances on Nd 143 is relatively negligible because, 
from the approximate Z{Z — 1)A~ 4//3 dependence expected from the Bethe- 
Weizsacker formula, 8E\ ^,SE^ 49 ~ 0.1 eV which is small compared to the 
near-threshold resonances in Nd 143 , specifically the one below the threshold 
at £ r 143 ~ -6 eV.] 

In deriving the bound (f49[) , we have implicitly assumed that, during the 
Oklo phenomenon, a took some fixed value a° kl ° (possibly different from 
a now ). The situation would be more complicated if, at the time, a(t) were 
oscillating on a time scale smaller than the duration of Oklo. As Sm 149 data 
depend essentially on the value of cr 14 9 at the end of the fission reaction, we 
expect that the bound ([49l) restricts the amplitude of the variation of a in 
many extended scenarios comprising a-oscillations. 

Dividing ( f49"|) by the age of the Oklo phenomenon, we can convert it into 
a bound on the time derivative of a averaged over the time span separating 
us from the end of the Oklo phenomenon. [In scenarios where a varies on 
the Hubble time scale, this averaged time derivative is nearly equal to the 
present time derivative.] This conversion introduces a further uncertainty, 
because the age of Oklo is not determined with precision. The geochrono- 
logical studies suggest an age around 1.8 x 10 9 yr |TB|, while several studies 
based on nuclear decay time scales gave ~ 10% higher values: for instance, 
1.98 x 10 9 yr g|, 1.93 x 10 9 yr gTJ, and 2.05 x 10 9 yr @. To remain conserva- 
tive in our bounds, we shall use the lower, geochronological value. Dividing 
Eq. (^) by 1.8 x 10 9 yr, we get the following conservative (95% C.L.) limit 
on the time derivative of a averaged over the time since the Oklo reactor was 
running 

- 6.7 x lO^V" 1 < - < 5.0 x lO^V -1 . (50) 
a 

This is weaker than Shlyakhter's estimates (which ranged between ±5 x 
10 -i8 yr -i m and iiQ-^yr-i |]| |TJ) but rests on a firmer experimental 



basis. On the other hand, this is between two and three orders of mag- 
nitude stronger than the other constraints on the variability of a (see the 
Introduction). Thanks to recent advances in atomic clock technology, it is 
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conceivable (and desirable) that direct laboratory tests might soon compete 
with the Oklo bound 

We have focussed in this paper on the time variation of the fine-structure 
constant because one can estimate with some confidence the effect of a change 
of a on resonance energies. It is more difficult to estimate the effect of a 
change in the Fermi coupling constant Gp, or, better, in the dimensionless 
quantity (3 = G F m 2 v cjh i ~ 1.03 x 10~ 5 . The estimates of Ref. @ for 
the (Weinberg-Salam) weak-interaction contribution to nuclear ground state 
energies yield -E^eak — ^weak — 5.6 eV. If one assumes that this gives an 
approximate estimate of the difference involving the relevant excited state of 
Sm 150 , and that there is no cancellation between the effects of changes in a 
and (3, one finds from Eq. the approximate bound 



I oOklo Onow I 

l - j^<0.02, (51) 

< KT 1 ^ 1 . (52) 



This bound is more stringent than the limit |/3//3| < 10 10 yr 1 obtained 
from the constancy of the K m decay rate @, and is comparable to the limit 
derived from Big Bang nucleosynthesis: \/3 BBN - f3 now \/(3 < 0.06 (H . 



Deducing from Oklo data a limit on the time variation of the "strength 
of the nuclear interaction" poses a greater challenge. First, one must notice 
that, as remarked in Section 2, only dimensionless ratios of nuclear quanti- 
ties, such as cr 149 /<7i43, can be extracted from Oklo data. Within the QCD 
framework, one generally expects any such dimensionless ratio to become a 
(truly constant) pure number in the "chiral" limit of massless quarks. 

Time variation of such a dimensionless ratio is then linked (in QCD) with 
possible changes in the subleading terms proportional to the mass ratios 
m q /m p , where m q denotes the masses of the light quarks. However, the 
chiral limit of nuclear binding energies is tricky because of non-analyticity 
effects in m q . The present chiral perturbation technology does not allow one 
to estimate the dependence of nuclear quantities such as El A9 or ai^/aus 
on m q /m p . One, however, anticipates that Oklo data might provide a very 
stringent test (probably at better than the 10 -7 level) on the time variation 
of m q /m p . To separate unambiguously the effects of variations in a and 
m q /m p , it would be necessary to extract from Oklo data several independent 



19 



measured quantities (e.g. by analyzing in detail the effects of resonance shifts 
in Gd 155 and Gd 157 ). 
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Table 1: Effective neutron cross sections of Sm 149 computed for 15 Oklo 
samples using published data. 
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Table 2: Oklo data corresponding to the extreme cross-section results of Eq. 
>0D. Notation as in Eqs. ©, (ITi 



Sample 
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Figure 1: Variation of the effective neutron capture cross section of Sm , 
0149, as a function of a possible shift A = E® Mo — E^ ow in the lowest reso- 
nance energy, for several values of the neutron temperature T. 0149, A and 
T (labelling the curves) are measured in kbarn, eV and degree Celsius, re- 
spectively. The two horizontal lines represent a conservative range of values 
of 0149 compatible with Oklo data. 
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